Optical Magnetoelectric Effect in Multiferroic Materials: 
Evidence for a Lorentz Force Acting on a Ray of Light 
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We theoretically propose that the optical analog of a Lorentz force acting on a ray of light is 
realized in multiferroic materials such as GaFeOa showing the magnetoelectric effect. The toroidal 
moment T = yV fj x Sj plays the role of a "vector potential" while its rotation corresponds to 
a "magnetic field" for photons. Hence, the light is subject to the Lorentz force when propagating 
through the domain wall region of the ferromagnetic or ferroelectric order. A realistic estimate on 
the magnitude of this effect is given. 
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The analogy between optics and particle dynamics has 
been a guiding principle in physics even leading to wave 
mechanics, i.e., quantum mechanics. In the classical the- 
ory, their analogy is established through the variational 
principle, which gives the Newtonian equation of motion 
for a particle and the eikonal equation of light. The 
variational principle for a classical particle is given by 
S J vols = with the velocity v = ^J2[E - U(x)]/m (E: 
total energy, U(x): potential energy, m: mass of the par- 
ticle). In geometrical optics, an equation for the trajec- 
tory of light is derived from the Fermat's principle; i.e., 
the time required to propagate from one to the other 
points is minimized. 5 J ds/V = S J dsn(r)/c = 0, where 
V = c/n is the velocity of light in the medium of re- 
fractive index n(f). This leads to the equation for the 
trajectory of light, which is similar to the equation of 
motion for a particle, as 



dj[ n (^w~] = Srad n(r). 



(1) 



In analogy to particle mechanics, we can interpret the 
above equation for light as follows. The term in the left- 
hand side is an "acceleration term" . On the right-hand 
side, the gradient of the refractive index Vn(f) acts as 
the "force" on the light and distorts the trajectory, giving 
rise to the refraction of light, i.e., Snell's law. Therefore 
n(r) is regarded as the "scalar potential" . In particle me- 
chanics, there is another type of force, i.e., the Lorentz 
force, in the presence of the magnetic field. This force 
is the transverse one to the velocity of the particle and, 
historically, has led to the concept of the "field" and rel- 
ativistic mechanics. Then the natural question to ask is 
whether the analogous force exists for light. It has been 
believed that there is no optical analog of magnetic field 
QQ, namely a Lorentz force on light in real space. On 
the other hand, a Lorentz force in fc-space was introduced 
recently in a polarization-dependent equation of motion 



as a generalization of the geometrical optics This 
"Lorentz force" in momentum space is described by a 
Berry phase in analogy to an anomalous Hall effect Q 
or a Magnus effect 4]. Haldane and Raghu focus on an 
edge state in a magnetic photonic crystal on the anal- 
ogy of the quantum Hall effect The edge state is a 
polarization-dependent one-way waveguide with broken 
time-reversal symmetry. However the Lorentz force on 
light in real space has been missing up to now. 

In this Letter, we theoretically propose that the 
Lorentz force on light independent of its polarization is 
realized in the multiferroic materials showing the non- 
reciprocal optical effect. The polarization-independent 
non-reciprocal optical effect requires a breaking of both 
inversion and time-reversal symmetries, which allows 
polarization-independent directional birefringence, called 
a magneto-chiral (MC) effect |6| |?J or an optical magneto- 
electric (OME) effect The difference between the 
MC and OME effects lies in their crystal structures — 
whether a magnet has a helical structure or a polar struc- 
ture. These effects have directional dependence, namely, 
k dependence in optical response. We focus on magnets 
with a polar structure, namely, multifcrroics. Novel op- 
tical properties in multifcrroics are characterized by a 
toroidal moment T = f% x Si 1 1 0l ] , where fi and S% are 
the displacement of the center position from atoms and a 
spin of an electron at ith site, respectively. We note that 
the definition of the toroidal moment implies that even 
in a antiferroelectric and antiferromagnetic system, the 
toroidal moment can be non zero. In ferroelectric and fer- 
romagnetic systems, we can rewrite the toroidal moment 
as T oc P x M, where P and M are spontaneous electric 
polarization and magnetization, respectively. Hereafter 
we assume that multiferroics are simultaneously ferro- 
electric and ferromagnetic for simplicity. We can reverse 
the toroidal moment by an external field by reversing the 
electric polarization or the magnetization. The toroidal 
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moment brings about the OME effect that refractive in- 
dex depends on whether the propagation direction is par- 
allel or antiparallel to the toroidal moment. The OME 
effect is expressed as — oc k ■ T, where arrows in 

the subscripts represent directions of light and k = . 

I 

We demonstrate below a novel refraction phenomenon 
caused by an inhomogeneous toroidal moment in mul- 
tiferroics. The equation of motion for light in multi- 
ferroics contains a Lorentz force in real space, where a 
spatially inhomogeneous toroidal moment gives rise to a 
magnetic field for light. Such spatial modulation of the 
toroidal moment can be achieved in a domain wall in mul- 
tiferroics. From optical viewpoints, the toroidal moment 
plays the role of a 'vector potential of light. 

Since the OME effect is independent of polarization, 
we can use geometrical optics, which is a scalar represen- 
tation, to describe the wave propagation through multi- 
ferroics. The OME effect is described by 



n(r) = no(r) + ak ■ T(r), 



(2) 



where k — k/\k\ is now represented as k = 4|. Substi- 
tuting the refractive index into Eq. (1), we obtain the 
following non-linear equation, 



— [n (r)fc] + a-^ [(& • T(r))k] = grad n (r) 



ds 

+a(k ■ V)f (r) + ak x rotf(r). 



(3) 



Each term on the Eq. © can be interpreted as follows. 
The first term on the right-hand side represents a "poten- 
tial force" for photons which yields a conventional Snell's 
law. The second term is the differential of the toroidal 
moment along the propagation direction, whose physical 
meaning is explained later. Remarkably, the third term is 
nothing but a Lorentz force on light with a vector poten- 
tial T{f). The k dependence in the second and the third 
terms gives rise to non-reciprocal refraction, namely, one- 
way propagation. 

In realistic materials, the contribution of the Lorentz 
force term can be observed in a domain wall(DW) of the 
toroidal moment. In a toroidal domain- wall (TDW) re- 
gion, the toroidal moment varies as T — (0, Bx, 0), which 
is analogous to the vector potential of a uniform magnetic 
field B in the Landau gauge, as shown in Fig. Since 
the toroidal moment is roughly related to the electric 
polarization and the magnetization as T oc P x M, the 
toroidal moment can be reversed by the reversal of P or 
M by an external field. In the TDW, we have no(r) = no 
and the equation is explicitly written as 

hqx + aB — (yxx) — aBy = 0, 



ds 
n y 



aB—(xy 2 ) 
ds 



0. 



where the dot denotes a derivative with respect to s. 



We assume the initial conditions to be 



k = (cos 8, sm8, 0), and we can integrate the equa- 
tions; 

tiqx + aByxx — aBy — hq cos 9, 
noy + aBxy 2 = no sin#. 

Then the solution to these equations is 



Ba 



2uq sin 



-x + - 



(4) 



s=0 



(0,0,0), 



which is interpreted as follows. The Lorentz force term 
of the right-hand side on Eq. J3J is explicitly writ- 
ten to be A: x rot T = (k y B, — k x B, 0), which gives 
rise to a cyclotron motion. The second term gives 
(k ■ V)T = (0, +k x B, 0) which cancels the y compo- 
nent of the Lorentz term. Therefore the right-hand side 
in Eq. ^ is proportional to (k y B, 0,0), leading to the 
parabolic solution. When the direction of the initial ve- 
locity is reversed, a sign of a curvature of the parabola 
is changed, which means that we have a one-way trajec- 
tory through a TDW. Because the toroidal moment is a 
physical field for optical waves and is related to the elec- 
tric polarization and the magnetization, it has no gauge 
invariance and can be controlled by external fields or tem- 
perature. It means that one can design a vector potential 
at one's disposal by spatially modulating P or M. 

Let us estimate the displacement 5 in Fig. 1. We put 
the thickness and the height of a sample to be I and h, 
respectively. The displacement is roughly estimated as 
5 ~ (h 2 /l)a, where a is the strength of the OME effect; 
a = {n^ — n t _)/no. In GaFeOa, these parameters are 
I ~ 100 nm, h ~ 10 fim and a ~ 10 -3 . Using these pa- 
rameters, we obtain the displacement to be 5 ~ lmn|llj. 
which can be observed through a microscope. A condi- 
tion for large S is to have a large amount of magnetic field 
and a large medium. Such a strong magnetic field can be 
achieved in a material with a large OME effect and a thin 
TDW, namely large a and small /, respectively. A large 
OME effect gives rise to a large refraction, and small I 
corresponds to a large differential of the toroidal moment, 
resulting in a large displacement. Here the validity of the 
geometrical optics should be discussed. Since the thick- 
ness of the DW in GaFeOa is lOOnm which is shorter than 
the wavelength of the visible light, the refraction of the 
light would be obscure. Therefore we should use ultravi- 
olet beam. However, a material with a DW whose thick- 
ness is larger than the wavelength of the incident light 
might be a good candidate to observe this non-reciprocal 
refraction. 

As an example, we consider an experimental set up 
for the TDWs shown in Fig. [2J Here, we introduce the 
TDW as a magnetic DW with uniform electric polar- 
ization. The incident optical ray is refracted through 
each DW, leading to two distinct regions as contrasts in 
a screen. One is a dark region represented by dark blue 
where photons cannot reach. The other is a bright region 
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characterized by a white color where both refracted light 
and transmitted light can reach. It means that TDWs 
can be observed as contrasts in a screen The re- 

gion in which the light propagates out of TDWs is repre- 
sented by a bright blue color. The contrasts purely have 
their origins in the spatially varying toroidal moment. 
To examine whether the contrasts are really caused by 
the inhomogeneous toroidal moment or not, we propose 
following two experimental methods. One is to introduce 
the light from another side of the sample. There should 
be different contrasts from those of another side, repre- 
sented by two screens 1 and 2 in Fig. [21 The other is to 
apply the external magnetic field ~ 500 Oe parallel to z 
direction to order the magnetization. The magnetic DW 
and the TDW vanish, resulting in the disappearance of 
the contrasts. Both of these two experimental tests are 
needed to confirm our proposal. 

The OME effect is usually very small for direct ob- 
servation: n_» - ~ 10 -3 for GaFe0 3 8), ~ 10~ 6 for 
Cr 2 3 13], - 10~ 4 for Eri. 5 Yi. 5 Al 5 0i2 .14] and so on. 
To overcome the disadvantage, one can enhance the OME 
effect in artificial structures. In periodic structures com- 
posed of multiferroics, the OME effect can be magnified 
by a factor of 1000 in multiferroic gratings and photonic 
crystals 0]. In contrast, what we demonstrate here pro- 
vides another example for a sizable OME effect. The 
non-reciprocal refraction obtained in this Letter is of the 
order of 1 /im, which can be observed without any use 
of a cavity, a photonic crystal or a patterned structure. 
Therefore, non-reciprocal refraction due to a toroidal mo- 
ment modulation provides us with a new method to mea- 
sure the OME effect. 

Kida et. al. found an OME effect in a submicron pat- 
terned magnet in which inversion symmetry is artificially 
broken by a chevron-shaped structure . The signal of 
the OME effect in the patterned magnet is of the or- 
der of 10 -3 , which is significantly large. They show that 
we can introduce a toroidal moment in a magnet with 
an artificially asymmetric structure. It may be possible 
that an inhomogeneous toroidal moment introduced by 
an artificial structure yields a Lorentz force on light. 

In conclusion, we have investigated a novel non- 
reciprocal refraction in multiferroics, which is polariza- 
tion independent. The toroidal moment, inherent in mul- 
tiferroics, plays the role of vector potential for light, and 
its spatial variation gives optical analogue of a Lorentz 
force. Photons feel no Lorentz force in a conventional 
sense under a real external magnetic field because they 
have no electric charge. The toroidal moment is a phys- 
ical field without gauge degrees of freedom, and can be 
experimentally controlled by an external field or temper- 
ature. It means that the Lorentz force can also be con- 
trolled by them. Such an interesting material as a mul- 
tiferroic matter has the potential importance of a novel 
external field for photons; the Lorentz force obtained here 
is one example. This Lorentz force on light fills one of 



the most fundamental pieces which has been missing for 
photon-electron analogy. 
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FIG. 1: Non-reciprocal refraction through a toroidal domain 
wall represented by a yellow box. The green arrows represent 
toroidal moments. Different propagation direction has differ- 
ent optical path, represented by the lines with red arrows. 



M0 




t<»t 



■ } 

7 



screen 1 



±®±v 



®M } 



I I 



multiferroics 
V 

z t) *-X 

\ screen 2 



FIG. 2: Non-reciprocal refraction through some domain walls. 
Domain walls are observed as contrasts on a screen, repre- 
sented by deep blue, bright blue, and white regions. The 
two screens 1 and 2 with different contrasts correspond to the 
incident directions, the face and the back, respectively. 



* Electronic address: sawada@appi.t.u-tokyo. ac.jp 
' Electronic address: nagaosaOappi.t.u-tokyo. ac.jp 
[1] S. Datta, Electronic Transport in Mesoscopic Systems, 

(Cambridge University Press, 1995). 
[2] M. Born and E. Wolf, Principles of Optics, 7th edition 

(Cambridge University Press, 1999). 
[3] M. Onoda, S. Murakami and N. Nagaosa, Phys. Rev. 
Lett. 93, 083901 (2004). 



4 



[4] K. Yu. Bliokh and Yu. P. Bliokh, Phys. Rev. E 70, 026605 

(2004) ; JETP Lett. 79, 519 (2004}. 

[5] F. D. M. Haldane and S. Raghu, cond-mat/0503588 

[6] L. D. Barron and J. Vrbancich, Mol. Phys. 51, 715 
(1984). 

[7] G. L. J. A. Rikken and E. Raupach, Nature 390, 493 
(1997); G. L. J. A. Rikken and E. Raupach, Phys. Rev. 
E 58, 5081 (1998); C. Koerdt, G. Duchs and G. L. J. 
A. Rikken, Phys. Rev. Lett. 91. 073902 (2003); F. A. 
Pinheiro and B. A. van Tiggelen, J. Opt. Soc. Am. A 20, 
99(2003). 

[8] M. Kubota, T. Arima, Y. Kaneko, J. P. He, X. Z. Yu 
and Y. Tokura, Phys. Rev. Lett. 92. 137401 (2004); J. 
H. Jung, M. Matsubara, T. Arima, J. P. He, Y. Kaneko 
and Y. Tokura, Phys. Rev. Lett. 93 037403 (2004); T. 
Arima, J. H. Jung, M. Matsubara, M. Kubota, J. P. He, 
Y. Kaneko and Y. Tokura, J. Phys. Soc. Japan 74, 1419 

(2005) . 

[9] S. Murakami, R. Shindou, N. Nagaosa and A. S. 
Mishchenko, Phys. Rev. B 66, 184405 (2002). 
[10] Yu. F. Popov, A. M. Kadomtseva, G. P. Vorob'ev, V. 
A. Timofeeva, D. M. Ustinin, A. K. Zvezdin and M. M. 
Tegeranchi, J. Exp. Theor. Phys. 87, 146 (1998). 



[11] For GaFeOa, optical path in Fig. 1 is not the case. The 
thickness / is so small in GaFeOa that optical beam is 
refracted and out of a TDW. The trajectory of the beam 
is a parabola in a TDW and becomes a line out of a TDW. 

,2 

The naive estimation S ~ — a breaks down. However, 
careful calculation suggests that the displacement still 
holds 8 ~ 1/im. 

[12] These contrasts simply represent whether photons reach 
or not, regardless of interference due to difference in light 
paths of refracted one and unrefracted one. Coherent 
light may make the interference pattern on a screen in 
addition to the contrasts. 

[13] B. B. Krichevtsov, V. V. Pavlov, R. V. Pisarev, and V. 
N. Gridnev, J. Phys. Condens. Matter 5, 8233 (1993). 

[14] G. L. J. A. Rikken, C. Strohm and P. Wyder, Phys. Rev. 
Lett. 89, 133005 (2002). 

[15] K. Sawada and N. Nagaosa, Appl. Phys. Lett. 87, 042503 
(2005). 

[16] N. Kida, T. Yamada, M. Konoto, Y. Okimoto, T. Arima, 
K. Koike, H. Akoh and Y. Tokura, Phys. Rev. Lett. 94, 
077205 (2005). 



